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Most of the studies on optimal transport are done for steady state regime conditions. Yet, there 
exists numerous examples in living systems where supply tree networks have to deliver products in 
a limited time due to the pulsatile character of the flow. This is the case for mammals respiration 
for which air has to reach the gas exchange units before the start of expiration. We report here 
that introducing a systematic branching asymmetry allows to reduce the average delivery time of 
the products. It simultaneously increases its robustness against the unevitable variability of sizes 
related to morphogenesis. We then apply this approach to the human tracheobronchial tree. We 
show that in this case all extremities are supplied with fresh air, provided that the asymmetry is 
smaller than a critical threshold which happens to fit with the asymmetry measured in the human 
lung. This could indicate that the structure is adjusted at the maximum asymmetry level that 
allows to feed all terminal units with fresh air. 

PACS numbers: 87.18.Wx , 89.75.Da , 89.75.Fb , 89.75.Hc 
Keywords: tree, asymmetry, optimal transport, lung, pulsatile flow 



Branched transportation networks are ubiquitous in 
living systems. Such trees have been suggested to allow 
efficient feeding the body volume from a small source, the 
aorta for blood or the mouth for oxygen [T]-[7] . For exam- 
ple, the human tracheobronchial tree brings fresh air to 
the oxygen-blood exchange units, called the acini, that 
fill the majority of the volume of the thoracic cage [1]. 

In the general perspective of how statistical physics 
may help to better understand the relation between 
structure and physiological function, we focus here on 
pulsatile trees in which the delivery of products has to 
be achieved in a limited time. For example in mammals 
respiration, the respiratory cycle is made on two succes- 
sive steps, inspiration and expiration. The period of this 
cycle is about 5 seconds for humans at rest (2 s for inspi- 
ration, 3 s for expiration). And of course, the transit time 
from the mouth to the acini has to be short enough so 
that expiration does not start before the arrival of fresh 
air into the acini. The present work analyses the possible 
statistical constraints related to this last condition. 

From the point of view of ventilation, the conduct- 
ing airway system can be modelled as an arrangement of 
pipes defined by their diameter and length. The branch- 
ings are essentially dichotomous, each airway being di- 
vided into two smaller daughter airways. A branching 
defines the beginning of a new generation. The tracheo- 
bronchial tree starts at the trachea (generation 0) and 
ends in the terminal bronchioles (around generation 15) 
at the entrance of the acini [l] . Since no gas exchanges 
take place in the tracheobronchial tree, its volume is re- 
ferred to as the dead space volume (DSV). Its value is 
around 170 mL in the human lung [8j. To better under- 
stand human respiration in relation with the tree struc- 



ture, we discuss successively the general properties of two 
models, respectively symmetric and asymmetric, before 
discussing the real human airway system. 

The first model has symmetrical branching. Although 
the human airway system is both non symmetric and 
exhibits some type of randomness, a basic step in de- 
scribing its morphology has been the introduction of the 
so called Weibel's "A" model [1 . In this first model, 
the tree is likened to a hierarchical network of cylindrical 
pipes with symmetrical branching and a uniform scal- 
ing ratio ho = 2~^/^ ^ 0.79 between the airway sizes of 
consecutive generations. This value corresponds to the 
classical Murray-Hess law |10] for which the diame- 
ter of the mother branch do and the diameters of both 
daughter branches di and d2 are linked by the relation- 
ship do = dl ^ d^. 

An important second step in modeling the human lung 
morphology has been the introduction from anatomical 
studies of a systematic asymmetry pTHT4] . This is the 
second model studied in this paper. The branching asym- 
metry is characterized by two different scaling ratios, 
ho,max = 0-88 and ho^min = 0-68 [14]. Each parent air- 
way gives rise to a larger daughter airway (the major air- 
way) and a smaller daughter airway (the minor airway). 
Note that (0.88)^ + (0.68)^ = 1 so that the asymmetric 
tree respects the above dissipation requirement. And so, 
the Weibel's A model is the exact symmetrisation of the 
real asymmetric structure. 

Those are the two models to be compared, includ- 
ing the possible role of a statistical noise associated to 
anatomical variability. Note that two different types of 
structural randomness can appear in these models: the 
first source of disorder, due to the systematic branching 



asymmetry, is found in the random succession of large 
and small airways along any given airway path. The sec- 
ond type of disorder (the variability) results from the 
statistical noise of the biological growth process [8 and 
spreads the distribution of the branch sizes at a given 
generation even in the symmetric model. 

The criterion that will be used in the following to quan- 
tify the delivery performance of the tree is the distribu- 
tion of oxygenation times of fresh air into the acini. The 
oxygenation time in one acinus is obtained by substract- 
ing from the total duration of the inspiratory phase, tins^ 
the time spent in the extrathoracic airways, text (approx- 
imately constant and equal to 0.47 s at rest [IS ), and the 
transit time from the trachea to this acinus, ttr- 



The performances of symmetric and of asymmetric 
trees have been computed on 15 generations trees. The 
geometrical parameters used in our computations, are 
summarized in Table [T[ These are the scaling ratios at 
each generation and the length to diameter ratio for each 
generation. It has to be noted that the values of the dead 
space volumes are kept almost similar. Due to the uni- 
form motion of the thoracic cage, each acinus is assumed 
to act as an hydrodynamic pump draining the same flux. 
In other terms, the gas exchange units are equitably ven- 
tilated [16]. As a consequence, starting from the bottom 
of the tree, any two daughters of a given mother branch 
create, independently of their sizes, the same additive 
flux in their mother branch. The time spent in a branch 
is then directly obtained from the flux and the branch 
size. In the inspiratory phase studied here, the total flow 
can be considered approximately constant in time with 
a velocity in the trachea of about 1 m/s [8 . Since the 
duration of the inspiration tins and the time spent in the 
extrathoracic airways text are the same for all acini, the 
oxygenation times in the acini are thus entirely deter- 
mined by the transit times of fresh air from the trachea 
to the terminal bronchioles. 

TABLE I: Model parameters 



Model 


Scaling ratio for D 


Ratio L/D ^ 


DSV (mL) 


Symmetric 


ho = 2-^rs 


3.00 


220 


Asymmetric 


ho,min = 0.68 
ho,max — 0.88 


3.00 


213 



^D and L: diameter and length of the airway. 

Fig. [l] shows the distribution of the acini oxygenation 
times for both models. The distribution for the symmet- 
ric tree trivially presents a single peak at t = 0.67 s as all 
the pathways from trachea to a terminal bronchiole are 
identical. The distribution of oxygenation times for the 
asymmetric tree is spread around a mean value of 0.82 s 
and has a standard deviation of 0.43 s. Therefore, on av- 
erage, fresh air arrives slighty sooner and remains longer 
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FIG. 1: Distribution of oxygenation times in the acini for 
both symmetric (red) and asymmetric (blue) trees. All the 
pathways have the same transit times in the symmetric case, 
hence a single peak distribution of oxygenation times at 0.67 s. 
In the asymmetric case the distribution is spread with a mean 
value about 0.82 s and a standard deviation of 0.43 s. 



in acini supplied by an asymmetric tracheobronchial tree 
than by a symmetric one. (Such an asymmetric tree 
structure can be described as multifractal [U]). 

We now study how the performances of these tree 
structures are robust or not with respect to anatomi- 
cal variability. To mimic this variability, we introduce a 
"growth noise" by adding Gaussian variations of the scal- 
ing ratios. At each branching, the values of the scaling 
ratios are modified to become random variables: 



hmin — ho^rnin (1 + O-X) 



(1) 



X being a centered gaussian random variable of standard 
deviation 1. The mean values ho^rnin and ho^rnax of these 
random variables correspond to the values given in ta- 
ble [l[ Scaling ratios of different bifurcations are assumed 
to be independent random variables while in the same 
branching, scaling ratios hmin and hmax are anticorre- 
lated. This means that if the random variable for hmax 
take a larger (resp. smaller) value than its mean value, 
then hmin is very likely to take a smaller (resp. larger) 
value than its mean value. 

The results are shown in Fig. |2] which displays the 
distributions of oxygenation times for both symmetric 
and asymmetric trees with randomized scaling ratios 
(cr = 5%). The distribution for the symmetric tree is now 
also spread, with an average oxygenation time of 0.67 s 
and a standard deviation of 0.28 s. Unlike the symmetric 
case, the distribution of times in the asymmetric model 
is not modified by the randomization of the scaling ratios 
(mean value of 0.82 s and standard deviation of 0.44 s). 
The delivery of products through a dichotomous tree of 
constant depth thus appears to be more efficient both for 
average oxygenation time and immunity versus growth 
fluctuations. 
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FIG. 2: Distribution of oxygenation times with randomized 
scaling ratios according to Eq. [l] symmetric tree in red and 
asymmetric tree in blue. The oxygenation times are spread 
in the symmetric case, but the distribution is shifted towards 
smaller oxygenation times. The distribution for the asymmet- 
ric tree is not modified by randomization. 



So far, the two model trees that we have considered 
had an equal number of generations for all pathways. In 
fact, the real human tracheobronchial tree is even more 
complex because not only the branchings are asymmetric 
but also the number of generations is not uniform [lT|[T2]. 
The human tracheobronchial tree does not end at a con- 
stant generation but at a constant airway diameter, that 
of the terminal bronchioles, around 0.5 mm \^16^. The 
terminal generations of the tree range from 8 to 22 [TT] 
and the length of the first generations airways also ex- 
hibits specific features linked to anatomical constraints 
[H [111 [18]. The diameter scaling ratios keep equal to 0.68 
and 0.88 for all generations. The aspect ratios (length 
over diameter) are specific for the first 4 generations (re- 
spectively equal to 3.07, 1.75, 1.43, and 1.85) and equal 
to 3.00 for higher generations (generations 6 to 22). The 
distribution of oxygenation times for the real human tree 
(not shown here) is found to have an average oxygenation 
time of 0.67 s and a standard deviation of 0.13 s. 

The asymmetry level, namely ho^rnax = 0-88 and 
ho,min = 0.68, used in the above computations was con- 
sidered as a given fact drawn from anatomical measure- 
ments. It is a natural question to ask for a reason of 
such values. For obvious reasons, the branching asym- 
metry cannot be too strong because it would lead to a 
structure with only very few wide pathways surrounded 
by a large number of much narrower pathways with large 
hydrodynamic resistances. Moreover, due to fiux conser- 
vation, the wider pathways correspond to larger transit 
times. If the asymmetry level were too important, a num- 
ber of extremities would not be supplied with fresh air 
because the transit time in their pathway would be too 
long. We will now investigate the influence of the asym- 
metry level on the acini oxygenation. The question that 
naturally arises is thus the following: how much asym- 



metry can there be? In other words, can one define an 
optimal asymmetry level? 

In order to investigate this question, trees of different 
asymmetry levels have been studied. In this study, the 
asymmetry level is characterized by one parameter a such 
that: 
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The scaling ratios 0.88 and 0.68 as measured in the hu- 
man tracheobronchial tree would for instance correspond 
to an asymmetry level of 36% {a = 0.36). All computed 
trees are conditionned to have the same thoracic volume 
(dead space volume + acini volume) and the aspects ra- 
tios L/D given above. Considering that fresh air has 
to remain at least 0.3 s in the acinar region in order 
to achieve the gas exchange process (a duration consis- 
tent with computations of the dynamical diffusion oxygen 
transport in the acinus [19 ), we have computed for each 
asymmetry level the proportion of acini with an oxygena- 
tion time larger than 0.3 s. This corresponds to a total 
transit time from the mouth to the terminal bronchiole 
smaller than 1.7 s. Results are presented in Fig. [3J The 
proportion of acini fed with fresh air for more than 0.3 s is 
100% for the symmetric tree (zero asymmetry level) and 
remains 100% until a threshold value of the asymmetry 
level. As one can see on Fig. [3J this threshold value is 
about 35%, almost identical to the value measured in the 
human tracheobronchial tree. 
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FIG. 3: Proportion of unactive acini (with oxygenation time 
smaller than 0.3 s) as a function of the asymmetry level ol. All 
acini are found to be active provided that the asymmetry level 
is below OL ^ 35%. This value almost exactly corresponds to 
the measured asymmetry level in the human lung (36%). 



So it seems that the systematic asymmetry found in the 
tracheobronchial tree corresponds to the maximum value 
that allows to feed all acini. The advantages of such a 
structure are several: the number of acini with fresh air 
is 100%, and the distribution of the ventilated volumes 
is robust against anatomical variability. Note that the 
spread distribution of arrival times in the exchange units 
may contribute to smoothen the oxygen delivery to the 



blood. On the other hand, due to this spread distribu- 
tion of transit times, ah acini do not receive the same 
volume of fresh air. Even if all acini receive the same 
flux, the uneveness of branching creates an inhomogene- 
ity of the volumes of supplied fresh air that multiplica- 
tively increases at each generation. In the limit of an 
infinite tree, this would mathematically lead to a multi- 
fractal distribution of the volume of fresh air delivered in 
the extremities [20]. This pre-multifractal behavior can 
be observed in Fig. [4] which shows a 2D representation of 
the distribution of volumes of fresh air delivered at gener- 
ation 10. One observes a wide spread of the distribution 
of fresh air volumes. This could indicate that a dynami- 
cal regulation of airway diameters might be necessary to 
minimize such effects, as already suggested to fight in- 
ertial effects in the upper part of the airway tree [21]. 
Note that the spread in the distribution of external gas 
in the acini would also induce an inherent noise in NMR 
imaging of the lung. 



"O CD 



CD O 




FIG. 4: 2D representation of the distribution of volumes of 
fresh air delivered at generation 10 in the human tracheo- 
bronchial tree (parameters described in the text). 

In summary, numerical computations show that sym- 
metric branching trees are not optimal to supply a vol- 
ume when a constraint of limited delivery time is im- 
posed. For trees of uniform depth, the average transit 
time is found to be smaller when the branching is asym- 
metric rather than symmetric. Moreover, the distribu- 
tion of oxygenation times in the acini of the asymmetric 
tree is almost not modified by a stochastic variability of 
the sizes of the branches. This approach is then extended 
to the study of the role of variable asymmetry in trees, 
with the same trachea, the same ratio length over diam- 
eter, the same diameter of the terminal bronchioles and 
the same inner volume. It is shown then that there ex- 
ists a maximum asymmetry level above which the num- 
ber of terminal units, or acini, supplied with fresh air, 
departs from 100%. Interestingly, this maximum value 
corresponds almost exactly to the asymmetry level mea- 
sured in the human lung. The geometry of the lung air- 



ways thus appears as if being adjusted to have the largest 
possible branching asymmetry, while still being able to 
feed efficiently all acini with fresh air. It may be consid- 
ered as remarkable that an asymmetric tree works better 
and that the natural selection of mammalian seems to 
have found a level of asymmetry that can be considered 
as best from the physical point of view. 

The authors would like to thank Pr E.R. Weibel for 
fruitful discussions. 
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